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Abstract: We present a subtraction scheme for eliminating the ultraviolet, soft, and collinear 
divergences in the numerical calculation of an arbitrary one-loop QCD amplitude with an 
arbitrary number of external legs. The subtractions consist of local counter terms in the 
space of the four-dimensional loop momentum. The ultraviolet subtraction terms reproduce 
MS renormalization. The key point in the method for the soft and collinear subtractions is 
that, although the subtraction terms are defined graph- by-graph and the matrix element is 
also calculated graph-by-graph, the sum over graphs of the integral of each the subtraction 
term can be evaluated analytically and provides the well known simple pole structure that 
arises from subtractions from real emission graphs, but with the opposite sign. 
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1. Introduction 

Tests of the Standard Model or one of its extensions frequently involve high momentum 
transfer processes in which one is looking for effects from new interactions or particles. In 
most cases, at least some of the particles interact via the strong interactions. Then, in 
order to obtain reasonably accurate predictions for the expected cross sections, it is necessary 
to calculate the cross section at next-to-leading order (NLO) in quantum chromodynamics 
(QCD). Sometimes this is possible using the currently available theoretical tools, sometimes 
not. This paper concerns a method for potentially extending the range of problems for which 
a next-to-leading order calculation is possible. 

An NLO calculation necessarily involves a virtual loop integration. That is, one has an 
integral 
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over a loop momentum / in d = 4 — 2e dimensions, with momenta ki, - ■ ■ ,kn leaving the 
graph at vertices around the loop. Of course, one wants d = 4 in the end, but there can be 
infrared divergences, which are temporarily regulated by working with d ^ 4. (There can 
also be ultraviolet divergences, but we ignore these for this initial discussion.) The infrared 
divergences arise because we work with a gauge theory containing massless particles. When 
some or all of the ki are lightlike, the integral for d = A can have divergences arising from 
regions of /-space in which the momentum of two of the propagators in the loop approaches 
a line parallel to one of the external momenta or in which the momentum of one of the 
propagators in the loop approaches zero. This leads to poles of the form l/e^ and 1/e in the 
integral. 

The traditional method of dealing with integrals like (|1.1| ), initiated in the context of 
collider physics cross sections by Ellis, Ross, and Terrano Q, is to calculate the integrals 
analytically. The result is expressed in the form of the residues of the and 1/e poles 
and a remaining finite piece, which contains the most important physical information. It 
would seem not to be helpful to have a result with and 1/e terms, but it is helpful. 
There are other integrations involved in the calculation of a physical observable. These are 
integrations over the momenta of final state quarks and gluons in Feynman graphs without 
virtual loops. As long as the observable has the property known as infrared safety, the 
integration over final state momenta will produce 1/e^ and 1/e poles that precisely cancel the 
poles from the virtual loop integrals. The integrations over final state momenta are generally 
much too complicated to perform analytically. One can subtract from the integrand certain 
simple integrands that match the complete integrand in one of the soft or collinear limits that 
gives a divergence. Then the difference is integrable in d = 4 dimensions and this integral is 
performed numerically. The integral of each subtraction term is added back, but this time 
with the integration performed analytically. This gives 1/e^ and 1/e poles that cancel against 
the poles from the virtual loop integrals. 

In order to make this method work, one needs analytical calculations of the virtual loop 
integrals that appear. Ellis, Ross, and Terrano [jjl supplied the one-loop graphs with n = 2, 
n = 3, and n = A external legs that occur in a calculation for e"*" + e~ — > 3 jets. Ellis 
and Sexton ||2| supplied the additional n = 2, n = 3, and n = 4 graphs that are needed 
for p + p ^ 2 jets. More recently, there have been heroic calculations by Bern, Dixon, and 
Kosower ^, |5|, ^ and by Kunszt, Signer, and Trocsayni for n = 5. With these results, 
one can calculate, for instance, +e~ — > 4 jets and p+p 3 jets. There are also calculations 
available with massive particles, such as massive quarks. 

One wonders if we could not get beyond n = 5 by simply performing the virtual loop 
integrals numerically. Not only might one be able to get to diagrams with more external 
legs, but also one could have the flexibility to change the lagrangian, add masses, make 
approximations, and so forth by simply altering the integrand. Of course, one is not going 
to perform integrals in d = 4 — 2e dimensions numerically. Instead, one would have to do 
something about the infrared divergences first, then take d = 4. 

Something like this was tried in Ref. |, g |lO|, |lll in the case of e"*" + e ^3 jets. 
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There, the approach was to perform the integral over the energy in Eq. (1.1) analyticahy 
and then to put all of the integrals over loop momenta / and final state particle momenta 
p together, summing over different contributions to the observable inside the integration. 
Then the real-virtual cancellations alluded to above happen inside the integrals, so that the 
integrals are convergent in three space dimensions. All of the integrals can be performed at 
once by numerical integration. 

In this paper, we present a different proposal, one that is very close to the traditional 
method. One would subtract from the integrand ( |1.1D certain simple integrands that match 
the complete integrand in one of the soft or collinear limits that gives a divergence. Then 
the difference would be integrable in d = 4 dimensions and this integral would be performed 
numerically. The integral of each subtraction term would be added back, but this time with 
the integration performed analytically. This gives the same 1/e^ and 1/e poles that one has 
in an analytic calculation, ready for cancellation as usual. 

We should mention that performing the 4-dimensional loop integral numerically is not 
completely without difficulties. In the massless theory, the integrand contains singularities 
1 / (p^ + iO) that lie on cones in the loop momentum space. One must deform the integration 
contour to avoid these singularities, as indicated by the "+iO" prescription. However, this 
same problem occurs in Refs. M, 0, |l^, 11 1, where it is handled in a straightforward manner. 



The issue addressed in this paper is to find what one should subtract from the integrand. 
We treat this problem in QCD with massless quarks, leaving the problem of other lagrangians 
and non-zero masses for future work. We find that there is a set of subtractions that is quite 
simple in the sense that the subtraction terms can be easily generated by the same sort of 
computer algebra that generates the original integrand in (LI). Furthermore, adding back the 
subtraction terms is straightforward: the residues of the 1/e^ and 1/e poles are the familiar 
ones, while the finite term that comes along with the poles is simple. 

There has been other recent work aimed at automating the calculation of loop integrals 
in such a way that those integrals that are not amenable to analytic evaluation could be 
evaluated numerically 13, 14 1. This approach makes use of the Feynman parameter 
representation of the diagrams and is quite ambitious in that it applies to multiloop diagrams. 
In the present paper we deal only with one-loop diagrams, although with an arbitrary number 
of external legs. Our treatment is then to subtract from the integrands certain terms that 
remove divergences, but to otherwise leave the integrand entirely as given by the momentum 
space Feynman rules. Thus we aim at maximal simplicity of the method. 



2. Strategy 

We seek a method for performing subtractions on one-loop graphs in such a way that we are 
left with an integral that can be performed numerically in four dimensions. That means that 
we must take care of ultraviolet, collinear, and soft divergences. 

We note first of all that in the case of tree graphs it is common to sum all of the graphs 
that contribute to a given set of external particles. Many graphs may contribute, but one 
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can sum these graphs algebraicahy before entering into a numerical computation. With loop 
graphs, the situation is different. If we are going to integrate over the loop momentum 
numerically, we (or our computers) have to know the structure of the integrand, which is very 
different for different graphs. Thus we will work graph by graph. 

Working graph by graph, the choice of gauge matters. We take the simplest choice, the 
Feynman gauge. The external lines in each graph represent physical, on-shell particles. Thus 
our algorithms are based on having spinors u{p,s), u{p,s), v{p,s) and v{p,s) for external 
quarks and antiquarks and polarization vectors £^{p, s) for external gluons. The polarization 
vectors can be defined using any physical gauge one cares to use. Our analysis makes use of 
the relation p ■ e(p, s) = 0. If it were desired, one could use unphysical gluon polarizations, 
with s)£'^(p, s)* — > —g^^- Then the relation p ■ £{p,s) = fails, and the reader is 

asked to substitute the following argument. We calculate a one-loop amplitude, which must 
be multiplied by a tree level amplitude to calculate a term in a cross section. Where p^ 
appeared contracted with a polarization vector, it now appears contracted with the tree level 
amplitude. But the contraction of p^ with the tree level amplitude vanishes. 

Thus our aim is to examine one- loop graphs with on-shell massless external particles. For 
each such graph, there may be ultraviolet, soft, and collinear divergences. In Sees. ^ ^, and ^, 
respectively, we will define subtractions for the integrand that eliminate each divergence. The 
ultraviolet subtractions are required for modified minimal subtraction (MS) renormalization. 
All that we have to do is to express the MS counterterm as a subtraction on the integrand. 
In the case of the soft and collinear subtractions, we do not want to change the value of 
the graph, so we add back what we subtracted. However, we subtract something from the 
integrand and add back its integral summed over graphs. What we need to add back is not 
simple graph by graph. However, when we sum over graphs we get a simple prescription for 
what to add back. We begin our investigations in the next section by setting up some useful 
notation. 

3. Matrix elements at tree and one-loop level 

In this section, we set up our notation. First, we define a notation for the color and spin 
structure of matrix elements, following in general the notation of Catani and Seymour jl5| ]. 
Then we add some notation specific to one-loop graphs. Finally, we define the wave function 
renormalization factors needed to relate amputated Green functions to the 5-matrix. 

The matrix element with m external QCD partons has the following general structure: 

M^''-''^''^'-''^{p^,...,Pm), (3.1) 

where {ci, . . . , Cm}, {si, . . . , Sm} and {pi, . . . ,Pm} are respectively color indices (which take 
— 1 values for each gluon and Nc values for each quark or anti-quark), spin indices (which 
take d—2 values for gluons in d-dimensional spacetime, and 2 values for quarks and antiquarks) 
and momenta. We take the momenta pi to be coming out of the graph, so that pf = with 
p^ > for outgoing partons andp-* < for incoming partons. The matrix element is computed 
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by contracting amputated Green functions with u{p,s) for an outgoing quark, v{p,s) for an 
outgoing antiquark, and e{p, s)* for an outgoing gluon. For incoming partons, the construction 
is analogous, using u{—p,s),v{—p,s) and e{—p,s). (There are identities available |l^ to 
simplify the spin wave functions, but these do not concern us in this paper.) 

It is useful to represent the matrix element as a vector in color spin space. One 
introduces a basis { |ci, . . . , Cm) ® l^i, . . . , s^)} in color (g) spin space in such a way that 

M'=i'-''=™;'*i'-'"™(pi,...,p„) = (^{ci,...,Cm\ {si,...,s^\^\Mipi,...,Pm)) . (3.2) 

With this notation, the matrix element squared (summed over final-state colors and spins) 
|Mp can be written as 

|M|2 = {M (pi, . . . ,Pm)\ M {pu ■ ■ ■ ,Pm)) . (3.3) 

3.1 Color space 

Following the notation of we associate a color charge operator Tj with the emission of a 
gluon from each parton i. Its action onto the color space is defined by 

T, = r,'',|a>|4/6| , (3.4) 

where is the appropriate color matrix {if cab if the emitting particle z is a gluon, if the 
emitting particle i is a quark and t^^ = —t^^ if the emitting particle i is an antiquark). That 
is, 

{ci, . . . , Ci, . . . , Cm, c\Ti\bi, . . . ,bi, . . . ,bm) = Sc^bi ■■■Tc^b,---^cmbm ■ (3-5) 

One also defines the corresponding adjoint operators t| for absorbing a gluon onto parton i. 
Then one defines a dot product of operators, Tj • Tj = T|Tj for emitting a gluon from line i 
and absorbing it on line j. For i ^ j this is 

T, • T, = T, • T, = T,%T^^,^ \c,\ \c,)^ . (3.6) 

The case i = j is special: 

T,-T, = C, , (3.7) 

where Ci is the Casimir operator, that is, Q = Ca = Nc li i is a gluon and Ci = Cp = 
(N^ — l)/2Nc if z is a quark or anti-quark. 

Note that, by definition, each vector \M) is a color singlet. Therefore color conservation 
is simply 

m 

Y,T^i\M{pi,...,Pm)) = ^. (3.8) 

i=l 

Using this notation, we also define the square of color correlated amplitudes, lA^^'^^p, as 
follows 

lA^^-'^P = {M{pi, . . . ,p™)| T, • Tfc \M{pi, . . . ,pm)) 

= 'M'^^-'-'>^---^{pu...,Pm)X n^aJb,a,M"------{pi,...,Pm). (3.9) 
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3.2 One-loop graphs 

As discussed in the introduction of this section, we wish to calculate the one-loop, m-parton 
matrix element graph by graph using numerical integration. In this subsection, we define 
some tools for this calculation. We use the vector notation in color and spin space of the 
previous subsection to write the one-loop matrix element as a sum of the one-loop Feynman 
graphs, 

\M{pi,...,Pm))= J2 \S{G;pi,...,Pm)) . (3.10) 

graphs G 

The 1^) vector is related to the one-loop, amputated Green function contracted with the 
appropriate spinors and gluon polarization vectors by 

gcu...,cm;su-,sm(^Q.p_^^^^^^p^-^ _ (^(ci,...,c„| (si,...,s„|)|g(G;pi,...,p™)> . (3.11) 

The graph G contains a one-loop, one-particle-irreducible subgraph T with n = n{G) 
external lines and n internal lines around the loop. Let us label momenta of the internal lines 
li,l2, ■ ■ - In- We integrate over the loop momentum. One can consider for some specific 
index k to be the independent loop momentum, which we call 1} Once we make a choice of 
Ik as the independent variable, any of the propagator momenta around the loop is related to 
/ by an equation of the form 

m 

k = lk + ^A{G,k),,p,, (3.12) 
i=i 

where the matrix A is determined by the structure of the graph and by our choice I = 1^. We 
write the integration over loop momentum as 




(3.13) 



Since the integration measure in the space of loop momenta does not depend on the choice 
of origin {e.g. dly = dlj.) we omit a label on dl indicating the definition of /. In the actual 
numerical integration, one chooses a definition for Z, then chooses a random value for / thus 
defined, and then computes the /j from Eq. ( |3.12| ). Thus our notation is 

\g{G-pi,...,pm)) = j^^\g{G;{l};pi,...,pm)) , (3.14) 

where Q denotes the integrand for graph G and {/} is a shorthand for /i, . . . , 

The loop integral could be ultraviolet divergent. This happens for n = 2, 3 and for 
n = 4 in the case that all of the lines entering the loop are gluons. In this case we define a 

^Alternatively, I could be the average of the Ik- This choice is convenient for constructing the ultraviolet 
counterterm for a loop that needs ultraviolet renormalization. See Sec. W. 
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suitable UV counterterm that reproduces the same result as MS renormalization. We define 
the counter terms in Sec. ^ 

In certain graphs, the loop integration is infrared divergent. For instance when an internal 
gluon line connects two external lines the integral has soft, coUinear, and soft-collinear poles. 
To handle the infrared divergences, we introduce simple counter terms that we subtract from 
the integrand for the numerical integration and then add back analytically. We discuss this 
in Sec. ^ for the soft divergences and in Sec. |6| for the collinear divergences. 

3.3 S-matrix 

The S-matrix can be constructed from the amputated Green functions by means of the LSZ 
reduction formula 

\Miull{pi,...,Prn))= (^V^^\M\2iPl,---,Prn)) • (3.15) 

Here we consider both sides of the equation to be vectors in color and spin space. The 
subscripts "full" indicate that this formula applies to the matrix element summed over orders 
of perturbation theory. The T (for "truncated" ) superscript on Ai on the right indicates that 
the matrix element is calculated by multiplying amputated Green functions by the appropriate 
Dirac spinors and polarization vectors. The factors \/r^, one for each external particle, are 
of two types, one for gluons and one for quarks and antiquarks. They are defined from the 
residues of the poles of the corresponding propagators at = 0. Specifically, when the fiavor 

the behavior 



of particle i is q or q then = R^, where the renormalized quark propagator near p"^ = has 



i5(p)-i?^^^ p^^O . (3.16) 
+ ilJ 

When the fiavor of particle i is g then = where the renormalized gluon propagator near 
p^ = has the behavior 

iD^^^ip) ~ Ra^^ + ■■■ p'^0 , (3.17) 

where the omitted terms are gauge terms proportional to p^ or p'^. 
The constants rj have perturbative expansions 

ri = l+ ^ + • • • , (3.18) 
where r^-^"* oc a^. Thus, including terms up to one-loop order, 

|-A^full(pi, . . . ,Pm)) = \Mtrcc{Pl, ■ ■ ■ ,Pm)) + \M^'^'^\pi, . . . ,Pm)) 

^rl^^\Mtree{pi,...,Pm)) • (3.19) 



-El 
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Here jA^C^'^)) is the one- loop contribution to Elsewhere in this paper we denote 

|_A4(r>i)^ simply as \M). 

Now we need to evaluate the residue constants r^^\ Consider first the quark propagator. 
Expanding to order Og we have 



iS{p) 



+ 



where is the one-loop renormalized quark self-energy graph. It has the form 

where evaluation of the graph in d — 4 — 2e dimensions gives 



Att e r(l - e) 



2 \<: 



r(i + 6)r(i-6)3(i-6) 

r(l - 2e) (1 - 2e) 



(3.20) 



(3.21) 



(3.22) 



and where Z^'^ is the one-loop MS counter term, 



.(1) 



^ 47r e r(l -e) 

Inserting this into Eq. ( |3.20| ), we find 



(3.23) 



(3.24) 



(1) 



We now want to take in order to find R^. The renormalization counter term 

removes the pole at e = for any negative value of . Having removed the pole, we can 
analytically continue e to Ree < 0. Then we have vl(p^) ^ as —p^ — > 0. This leaves 



iS{p) 



p^ + iO 



1-4"} ■ 



(3.25) 



Comparing to Eq. (3.16), we have 



R 



(1) 



-Z^^ = Cf ^ ^^^^ 



4^ e r(l-e) 
The analysis for gluons is similar. One has 



(3.26) 



(1) 



^(1) 



as 1 (47r) 



47r e r(l - e) 



(3.27) 



where Tr = 1/2 and ?^/ is the number of quark flavors. 
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4. Renormalization 



In the following sections, we will consider the divergences of infrared origin that can arise in 
the integral 



/ 



'^'^ \g{G;{l};pu...,Pm)) (4.1) 



(2vr) 



for a graph G. First, however, we should deal with the possible ultraviolet divergences. These 
are to be eliminated according to the MS renormalization prescription. However, since we 
calculate loop integrals by numerical integration, the implementation of the MS prescription 
needs some analysis. This is the subject of the present section. 

Consider an ultraviolet divergent one-loop graph with n propagators in the loop. Here n 
could be 2, 3, or 4. Denote the momenta leaving F by /ci, . . . , In our application, F will be 
a subgraph of G and the momenta ki, . . . ,kn will be linear combinations of pi, . . . ,pm- The 
graph F has the generic form 



F(fci, ...,kn)= I F(A:i, ...,kn;l) . (4.2) 



The functions F and F may carry spinor and vector indices as well as color indices. Our 
notation here suppresses these indices. We seek to calculate the renormalized version of F. 
With MS renormalization, this is 

[r]R = lim{F-[F]p„i,} , (4.3) 

where 

[rU = 7f|^ xlim[.F(.)] . (4.4) 

Here we are to choose the external momenta so that F does not have infrared divergences. 
For instance, all of the external momenta can be spacelike. Then the only pole present in 
F(e) in Eq. (4.4) is the ultraviolet pole that is to be removed by renormalization. 

Since we are performing integrals numerically, we need to represent [Fjp^^jg as an integral 

(2^« 



[r]poie= j:—iruv{ki,...,kn;l) (4.5) 



in such a way that can be calculated as 

^^-^ lim|r(A;i,...,A;„;/) -r;7V'(/ci,...,fc„;0} • (4-6) 

This last step is justified if the integrands F and Tuy match up to an l^^ remainder for I oo 
at fixed e and if Tuv if free of infrared singularities. 
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There is more than one possibihty for Tf/y. However, there is a simple prescription that 
works in all cases save one, the one-loop gluon self-energy, which has a quadratic divergence 
and needs a more elaborate treatment. The prescription, ignoring complications arising from 
the tensor structure of F, is as follows. (We give the full details including the gluon self-energy 
and the tensor structure in Appendix ^) . Write the starting integrand in the form 

[If + lO) • • • [li + lO) 

Here we display the denominator for each propagator and put everything else in a numerator 
function N{ki, . . . , k^, l). We make a definite choice for the loop momentum In general, 

li = l + Ki , (4.8) 

where Ki is a linear combination of the momenta kj leaving the graph at the vertices. We 
define 

l = -[h + h + --- + ln] ■ (4.9) 
n 

Then 

^K, = . (4.10) 

i 

The numerator contains an overall factor times a polynomial in /, with coefficients that 
can be polynomials in e. Except for the gluon two-point function, which we treat separately, 
the highest order term is of order corresponding to a linearly divergent integral for 

d = 4. In some cases, this term is absent. The next highest order term is of order 
corresponding to a logarithmically divergent integral. This term is present as long as we are 
dealing with a graph that needs renormalization. Call these two terms Njjy, 



N{ki, ...,kn;l) = Nuviki, ...,kn;l) + 0(/""-^) . (4.11 



Now define 



r^v(fei,...,fc.;0 = , (4.12) 

where [i is the MS renormalization scale and A is a constant to be determined. With this 
choice, Vuy is free of infrared singularities that could lead to an infrared divergence. It 
matches F for large I. Here it is evident that the leading terms, which correspond to the 
linear divergence, match. To treat the next terms, note that for I oo, 

■ ' ^'"■''■^O(^) . (4.13) 



il + Kiy---il + Kn)^ (P)" (/2)"+l '^\^(p)n+l 

The second term vanishes because of Eq. (lilol) . For this reason, the order ^ term in 
F — F[7y, which would give a logarithmic divergence, vanishes. Thus F — Tjjv falls off fast 
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enough that its integral is ultraviolet finite when d is near 4. Consequently, the e — > limit 
of the integral of F — Tuv can be evaluated by simply setting d = 4. 
We note next that the integral 

^Tuv{h,...,kn;l) (4.14) 

will be a constant times the MS pole factor (1/e) (47r)'^/r(l — e) plus a remaining finite piece. 
The constant multiplying the pole is necessarily the same constant as for the integral of F 
since the ultraviolet behaviors of F and F[/y match. The finite constant depends on A, so it 
is a simple matter to choose A such that the finite part vanishes, yielding Eq. ( ^.51 ). 

We have seen here how to construct the renormalization counterterm for a divergent 
subgraph F of a graph G. Using this counterterm, renormalization amounts to replacing 

^\g{G;{l};pu...,Pm)) (4.15) 

by 

{\G{G;{l};pi,...,p^))-\n{G;{l};pi,...,p^))} , (4.16) 

where \TZ{G; {l};pi, ■ ■ ■ ,Pm)) is obtained from G by substituting Fj/y for F. In the case that 
graph G is ultraviolet finite, we define \TZ{G; {l};pi, . . . ,pm)) to be zero. We will use this 
notation in subsequent sections. 

We give specific expressions for Tjjvi^i: ■ ■ ■ , ^n; I) in Appendix There, we also exhibit 
what to do with the tensor structure of the F(A:i, . . . , /) and with the gluon self-energy di- 
agram, which needs a somewhat different treatment than presented above because it contains 
a quadratic divergence. 

5. Soft singularities 

The integrand of a one-loop graph becomes singular when the momentum of an internal gluon 
loop line that connects to two external lines becomes soft. Consider a one- loop graph G with 
external momenta pi,P2, ■ ■ -Pm directed out of the graph. Choose two of the external lines, 
with labels i and j. If the lines i and j connect via three-point vertices to the two ends of a 
gluon propagator in the loop, then we will say that {i,j} is in the soft indices class for graph 
G, £ Is{G). In this case the integrand \Q{G; {l};pi, . . . ,pm)) is singular in the limit 

that the momentum of the gluon propagator that connects lines i and j, call it l^, tends to 
zero. The loop integration is logarithmically divergent at this point. 

The denominators that become singular in the soft limit 1^ ^ are for the soft gluon 
propagator, (Z^ +Pi)'^ for the immediately preceding propagator in the loop and (4 —VjY for 
the immediately following propagator^. It is useful to define a soft limit function 

\fij{G;Pl,. . . ,Pm)) = lim llih + Pifih - Vjf\G{G;ll, . . . ,ln;Pl, ■ ■ ■ ,Vrn)) ■ (5-1) 



^Here we choose the positive direction around the loop to be from line i to line j. 
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Here we implicitly take Ik to be the independent loop momentum and use Eq. (|3l^) to 
determine the other /j. With this definition, \fij) vanishes if ^ Is{G)- We can then 

define a soft subtraction for the integrand G as 

(/^ + lO)((/fc + Pi)^ + lO)((/fc - Pj)^ + lO) 

The integral for the original graph minus this subtraction is 

d'^l { ~ ~ 1 

\G{G\{l]\pi,. . . ,pjn)) -\Sij{G;lk;pi,. . . ,Pm))\ ■ (5.3) 



d 



(2vr) 

By construction, this integral does not have an infrared divergence from the region 1^ — > 0. 
It may have other divergences, which will be eliminated by other subtraction terms. For 
instance, there may be another pair of external lines that connect to the virtual loop 

and are joined by a virtual gluon line with label k' . Then we should also subtract as 
defined above, from the graph. There may also be collinear divergences, to be discussed in the 
following section. Finally, there may be ultraviolet divergences, as discussed in Sec. ^. The 
idea is that after all divergences are subtracted, the indicated loop integral can be performed 
by numerical integration. 

Having subtracted the integral of we should add it back. For this purpose, we need 
to study the structure of the soft subtraction in more detail. The first step is to extend the 
vector notation that we have been using to cover unphysical gluon polarizations. Previously 
we have implicitly supplied a projection operator |s)^,s| onto the d — 2 dimensional space 
of physical polarizations for each external gluon. For the purpose of the present section we 
consider a d-dimensional polarization space. 

In addition to the physical polarizations represented by |s) we introduce polarizations 
along the coordinate axes. These have the orthogonality relations 

(^k) = ~9t^y ■ (5-4) 

The minus sign is present because we have chosen the inner product in the bra-ket spin space 
so that {s'\s) = +5sis- The corresponding completeness relation for the |//) vectors is 

1 = -\^J)9^''{v\ (5.5) 

with an implicit summation. We take the sign convention for the vectors |^) such that 

{li\M) = . (5.6) 

Then our previous definition = e*{p,s)^M^ can be written 

{s\M) =e*{p,sY{^l\M) . (5.7) 

The soft limit function has a rather simple form 

|4(G;pi, . . . ,pm)) = -iglii^'Ap,-pjTi ■ Tj Hij\g{Gij{Gy,pi, . . . ,Pm)) , (5.8) 
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where the graph Gij{G) represents the tree level amputated graph obtained from graph G by 
omitting the three singular propagators and the vertices where they join the external lines i 
and j. Here Hij is an operator in the spin space of the partons i and j. The form of Hij 
depends on whether partons i and j are quarks (or equivalently antiquarks) or gluons. In 
computing Hij, we keep only terms that survive under projecting onto physical polarizations 
for the final state partons. We find 

Hij{qq) = 1 , 

m.igg) = 1 + 14 ^,(«| + kX- + \f.\\i^y^-^^ MM ■ (5-9) 

' * 2pi ■ Pj * ' ' 2pi ■ Pj ' I n\ /J ^p. .p. i\ \j\ I 

Notice that in each term beyond the unit operator there is a factor (p\ or pf .{^a|. This 
feature will play an important role in the subsequent analysis. 
With this result, the soft subtraction is 

\Sij{G;lk]Pi, . . . ,Pm)) = Eij{lk,Pi,Pj) Ti ■ Tj Hij \g{Gij{G);pi, . . . ,pm)) , (5.10) 



where the Eij{lk,Pi,Pj) is the familiar eikonal factor [|^, 19 1 



(ifc + iO)((/fc + iO)((/fc -PjT + 1O) 

Note that, in Eij, one could arrange the definitions so that the denominator (/^ +iO = 

2Zfc -pi + ll + iO is replaced by 21^ -pi + iO, as in [18|, since these expressions match in the limit 
Ik 0. Similarly (/^ — pj)^ + iO could have been —21^ ■ Pj + iO. The form used in Eq. ( ^.ll]) 
suits our purpose here because its integral over 1^ is ultraviolet convergent. The subtraction 
term in Eq. ( ^.lOD is illustrated in Fig. ||. 

We need to add back the sum over graphs of the integral \Sij) of \Sij). We are now 
prepared to calculate what this quantity is. Defining V|°^*(e) to be the integral of Eij{l), we 
have 



\Sij{G;pi,. . . ,pm)) = j j^^\Sij{G]{l},pi,...,pm)) 

= V!f{e) T, • T, Hij \G{Gi,{G)-pi, . . . ,Pm)) • (5.12) 
The integral can be performed analytically, with the result 

Kf'i^^^f£,E,i,-^^(^^(U0i.^ . (5.3, 



J {2'kY ^ AttTII- e) \-2pi-pj ) \e'^ ^ 

Thus \Sij{G)^ is a known integral times a tree level graph with a modified color factor and a 
kinematic factor Hij. 
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lim 





Figure 1: Pictorial representation of the soft gluon subtraction, Eq. (5.1C). The quantity depicted is 
{\l Ei.j)\Si.j'^ for i = q, j = g. The double dashed line represents the octet color exchange contained in 
the factor Ti-Tj ■ The two terms correspond to the two terms in Hij (qg) . In the second term, the arrow 
on the gluon line represents the factor pj ^.(/3|, which has the effect of replacing the polarization vector 
£*(j>j, s) that would normally be contracted with the amputated Green function by the momentum pj. 



The polarization vector is then contracted with pi 
in which {s\iy) Pi ^ -e*{pj,s) - p.. 



The minus sign arises from our sign conventions. 



The sum over graphs G of the subtraction terms is even simpler. We consider a given 
set of final state parton flavors and a given choice of The tree graphs Gij{G) have the 

same parton flavors and momenta. Summing over all graphs G for which {i,j} G /^(G), the 
graphs Gij{G) cover all tree graphs with this choice of external partons. Thus summing over 
graphs G gives the full tree amplitude |A^tree(pi; • • • ,Pm)) corresponding to the given final 
state. On the other hand the tree level matrix element is gauge invariant, that is 

^{n\Mtree{pi,---,Pm)) =0 , (5.14) 

where j is a gluon leg. Thus gauge invariance ensures that the complicated terms in the 
operator Hij, Eq. (|5.9| ), do not contribute, so that Hij can be replaced by the unit operator. 
Then the soft contribution that we need to add back is 

\S^J(.G■pl, . . . ,pm)) = V!f\e) Ti ■ \Mtree{pi, ■ ■ ■ ,Pm)) • (5.15) 

G 

We now summarize the soft gluon subtractions, this time including a sum over pairs of 
indices {i, j}. We subtract 

^ \Sij{G]lk;pi,. . . ,Pm)) , (5.16) 
{ij}e/s(G) 

defined in Eq. ( ^I^ ) from the integrand \Q{G; {l},pi, . . . ,pm)) for each one-loop graph G. 
Then we add the integrals of these terms back in the form 

V^f (e) T. • |A^tree(pi, ■ • • ,Pm)) ■ (5.17) 

{id} 

In Sec. 1^, we will see that the terms added back are cancelled by similar terms corresponding 
to real soft gluon emission. The terms subtracted from each one-loop graph serve to remove 
the soft gluon divergences. This leaves the collinear parton divergences, which are the subject 
of the next section. 
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6. Collinear singularities 



One-loop graphs have logarithmic infrared divergences that arise from integration regions 
in which the momentum on an internal loop line that connects to an external line becomes 
collinear with the momentum of the external line. In this section, we define a term that, when 
subtracted from the graph, eliminates the divergence. This construction is quite similar to 
what we needed for soft gluon divergences. However, the structure of the collinear subtraction 
term is more complicated, so we will need a more involved analysis. 

Consider a one- loop graph G with external momenta pi,p2, ■ ■ ■ Pm directed out of the 
graph. Choose one of the external lines, with label i. If the line i connects via a three-point 
vertex to the loop, and if the loop partons are both gluons or are one gluon and one quark 
or antiquark, then we will say that i is in the collinear index class for graph G, i £ Ic{G). In 
this case the loop integration has a logarithmic divergence arising from the region in which 
the momenta of the loop propagators that connect to line i, call them j and j + 1, become 
collinear to the outgoing momentum of line i, 

^ ^ Pi 

-Ij+i {l-x)pi , (6.1) 

with < X < 1. 

The denominators that become singular in the collinear limit (SJ.) are /| and (— /j+i)^ = 
{Pi — I'j)'^- The coefficient of x {pi — Ij)'^] in the integrand for the graph is non-singular in 
the collinear limit. It is useful to define a collinear coefficient function \ f^'^{G; x;pi, . . . ,Pm)) 

by 

|/f'°(G;x;pi,...,p„i))= lim ^ {pi - Ijf \Q{G;li, . . . ,ln;pi, ■ ■ ■ ,Pm)) ■ (6.2) 

ij- Pi 



In writing Eq. ( |6.7D , we implicitly take Ij to be the independent loop momentum and use 
Eq. ( |3.12| ) to determine the other li in \Q{G;li, . . . ,ln',Pi, ■ ■ ■ ,Pm))- In particular, /j+i = 
Ij - Pi- 
One can use the collinear coefficient function to remove the collinear divergence. Consider 
the integral 



(2^ 



\g{G;{l},pi,...,pm)) 



dx5[x-'-l-^]\ff^\G-x-p^,...,Pra)) \ . (6.3) 



,Pm)>| 



(/2 + i0)((p, -/j)2 + i0) Jo V Pi-rii 

Here we define the momentum fraction x away from the collinear limit by using a lightlike 
vector rii. A good choice for Ui is 



pf + ^^w^ , w^= Pk ■ (6-4) 



fcgfinal state 
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By construction, the integral in Eq. ( |6.3| ) does not have a divergence from the colhnear 
region ( |6.1|) . The divergence was only logarithmic and the subtraction removes the leading 
singularity, leaving at worst an integrable singularity. 

There are two problems with Eq. (|6.3[ ). The first is that we have already subtracted 
something from the integrand for the graph G, so we do not start with a clean slate. In the 
present notation, the soft subtraction term corresponding to propagator j in the loop being 
soft is, in the collinear limit (|6.1|), 



(/| + iO)((p. - l,f + iO) I i^-^o • • • '^-)> ■ (^-^^ 

(If the parton in loop propagator j is not a gluon, then this quantity is zero.) Similarly 
the soft subtraction term corresponding to propagator j + 1 in the loop being soft is, in the 
collinear limit (|6.1|), 

[i] + io)((p. - i,Y + io) i^^i " • • • '^-^^ • 

In order to cancel these contributions to the net integrand in the collinear limit, we should 
subtract them from the collinear subtraction term in Eq. (|6.3| ). To this end, we define a 
revised collinear coefficient function \ ff^ {G]X]pi^ . . . ,Pm)) 

\ff{G;x;pi,...,pm)) = . . . 

1 1 

lini y\f^'°{G;y;pi,...,pm)) - hm (1 - y)|/f '°(G; . . . . (6.7) 

With this coefficient function, our subtraction inside the loop integral is 

^ f\xs(x-^-^)\ff{G;x;pi,...,pm)) ■ (6.8) 



(Z|+i0)((pi-/,)2 + 

A second problem remains with Eq. (|6.8| ) . The integral of the subtraction term is ultra- 
violet divergent. We can easily fix that by modifying the subtraction term to be 



where 



fuv{h.h-P^) = \{Y, ^^ + 7^ ^.J • (6.10) 

Here ^ is the MS renormalization scale and e = 2.71828 ... is the base of natural logarithms. 
The factor fjjv provides an extra power of /| in the denominator for large Ij but equals 1 in 
the collinear limit. 
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In summary, we subtract from the integrand for each graph G colUnear subtraction terms 



J2 \CiiG;{l},pi,...,pm)) , (6.11) 



where 

\Ci{G;{l},pi,. . . ,pr, 



fuvi^j, - Pi) 



(/2+i0)((p,-/,)2+i0) 
X 



/■'dx^fx-^) |/f(G;x;pi,...,p^)) . (6.12) 
JO V Pi ■ nij 



These subtractions, together with the soft subtractions, remove all of the infrared divergences 
from the loop integrals. One should note that, although the definition, Eq. ( |6.7D , of the 
collinear subtraction terms seems rather complicated, it is a purely algebraic construction 
that can be accomplished by straightforward computer algebra. 

Our next task, pursued in the following two subsections, will be to add back the subtrac- 
tion terms \Ci{G; {l},pi, ■ ■ ■ ,Pm)), this time integrated and summed over graphs G. That is, 
we want to add back the sum over graphs of 

\C,{G-{p,,...,Pm))^ I ^ |C,(G;{/},pi,...,p^)) . (6.13) 

6.1 Quark line 

For each one-loop virtual graph G and for each of its external lines i with i G Ic{G), we have 
subtracted a quantity \Ci{G;pi, . . . ,Pm.)) from \G{G;pi, . . . ,Pm))- These subtractions soften 
the collinear divergences in \Q{G;pi, . . . ,pm)) with its soft gluon subtractions, so that one 
could perform the loop integral for graph G by numerical integration. Now we must add the 
collinear subtractions |Cj(G;pi, . . . ,Pm)) back, this time performing the integral analytically. 
By design, these integrals have collinear singularities that give poles at d = 4 dimensions. 
The idea is to cancel these poles against the poles from the subtraction terms for collinear 
divergences in NLO tree graphs. 

Unfortunately, the structure of the \Ci{G;pi, . . . ,Pm)) for any given graph G is quite com- 
plicated. Fortunately, the structure of the sum over graphs G of \Ci{G;pi, . . . ,Pm)) is vastly 
simpler. One gets a singular factor times the leading order matrix element |7W(j)i, . . . ,Pm)) 
calculated from tree graphs. To see this, we simply need to apply the construction used to 
prove factorization of collinear parton contributions in hadron-hadron collisions . 

In this subsection, we consider the case that line i is a quark line. The case that i is an 
antiquark line follows trivially and is covered at the end of this section. The case that i is a 
gluon line is treated in the following subsection. 

We begin by setting up a useful notation. By assumption, line i, carrying momentum pi 
out of the graph, connects to a virtual loop in G. Line i must therefore connect to the loop at 
a quark-gluon-quark vertex. We choose to label the propagators in the loop so that the gluon 
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Figure 2: Illustration of Eq. ( 6.15| ). The left hand side represents f^''^{G;x;pi 
defined in Eq. (6.2) as the collincar limit of {Ij 
Here the integrand G consists of a spinor u{pi 



,Pm)", which is 

Pi)2 times the integrand t/(G'; /i, . . . , Z„;pi, . . . ,Pm)°'- 
a vertex and two propagators times an amputated 



Green function V, as in Eq. (6.14). The illustration represents the integrand: an integration over the 
loop momentum is not implied. On the right hand side, the gluon line of V is contracted with (xpi)^, 
indicated by the arrowhead, while the quark line is contracted with u{(l — x)pi). The remaining factors 
in the right hand side of Eq. ( |6.15 ) are indicated. 



line in the loop has label j and the quark line has label j ' + 1. Thus a gluon line with label 
j carries momentum Ij into the vertex and the quark line with label j ' + 1 carries momentum 
—Ij+i = Pi — Ij into the vertex. We write the integral for the graph as 

f _ 

G{G,pi,...,PmT = I j:^giG-lu---,ln,Pl,...,Pmr 

d% 1 

(z|+io)((p.-/,)2+io) 

^u{pi)-1^{f.-tj)\y^'{Gi,pi,...,Pi-lj,...,Prn,lj)]l^, ■ (6.14) 

Here we have abandoned our vector notation for spin and color space and written the spin and 
color indices that are needed for the calculation explicitly. We display a quark color index a 
for quark line i m. leaving the other indices on Q unwritten. The amplitude V has a color 
index a corresponding to the gluon j in the loop and a quark color index a' corresponding 
to the quark line j + 1 in the loop. There is an explicit color matrix T^^, connecting the 
colors of V to the color of Q. The amplitude V carries a vector index jj, corresponding to the 
polarization of the gluon line j. It also carries a Dirac spinor index for the quark line j + 1. 
However, we use the matrix notation for the Dirac structure of this line without displaying 
the Dirac indices explicitly. The Dirac spinor u{pi) represents the final state quark on line i. 
We have displayed the integration over the loop momentum Ij , the quark-gluon-quark vertex 
and the propagators for the quark and the gluon in the loop. Everything else is included in 
the Feynman amplitude V for the amputated tree level graph Gi obtained by omitting the 
propagators j and j + 1 and the vertex that attaches these propagators to external line i in 
graph G. 

Thus the amplitude V is defined by Eq. ( |6.21| ) . The momentum arguments of V are the 
external momenta of the original graph, plus the momentum pi — Ij of the the quark line and 
the momentum Ij of the gluon line. 




With the notation thus defined, we are ready to calculate. The subtraction term f'-"'^ 
defined in Eq. ( |6.2D can be expressed in terms of V as 



f[''''^{G;x;pi, . . . ,pmT = lim I'j {pi - Ij)'^ g{G;li, . . . ,ln;Pi, ■ ■ ■ ,PmT 

Ij >X J) 2. 

= ^Ss^^^T^a'U{Pih,^{l - X)f. [V''(Gi,pi, . . . , (1 - X)pi, . . . ,Pm,XPi)]l, 



2ig,/iXa'^^^ {XP^)^U{{1-X)P^) 

x[V^'{Gi,pi,...,{l-x)pi,...,pm,xpi)]l, . (6.15) 



In the last step, we rewrite u{pi) as u{{l — x)pi)/\/l — x in order to have a spinor evaluated 
at the momentum of the corresponding quark line. This equation is illustrated in Fig. ^. 

We would now like to sum over graphs. We consider one-loop graphs G with m external 
particles having flavors {fi, . . . , fm} and momenta {pi, . . . ,Pm}- We are considering the 
collinear subtraction for parton i, which we assume here is a quark. Our graphs are amputated 
on their external legs. Furthermore, we need consider only graphs that have a collinear 
divergence for external leg i: i £ Ic{G). Let us call this class of graphs G. 

The corresponding graphs Gi are amputated tree graphs with m + 1 external legs. The 
flavors of the first m external particles are {/i, . . . , fm]-, the same as for the graph G. The 
momenta of these particles are the same as for G except for particle i, which carries momentum 
{\ — x)pi instead oipi. The external particle with index m + 1 is a gluon with momentum xpi. 
Let us call the set of all such graphs C . When we sum over all graphs G £ G, the graphs Gi 
include all graphs in G' except for the graphs in which gluon m + 1 couples directly to quark 
i. These graphs are absent because graphs with a self-energy insertion on external line i are 
not included in G. Let us call the set of (m + l)-particle graphs that we do get G'^ and the 
graphs that we don't get G'_- 

The collinear subtraction defined in Eq. ( |6.12| ), summed over graphs G G C is 

^C,iG;pi,...,pmr 

d'^lj fuv{lj,h-Pi) , J Ij-rii 



f dx5(x-^-i^)^f{G;x;pu...,Pmr 
Jo V Pi-niJ 



{27rY [1] + iO){{pi - IjY + iO) io K -n/ 

d% fuv{l,,l,-p^) rdx _ ^.^^^.^a^^ 



(27r)^ (/2 + iO)((pi - IjY + iO) Jq X V P^■n^, 

x{{l-x)H''^,{x-pi,...,Pm)-K'{^,p^,...,Pm) \ . (6.16) 



Here we define 



Ha'{x;pi, . . . ,Pm) = 

(l-x)-^/2 ^ q^u{{l-x)p,)[Y^'{Gi,pi,...,{l-x)pi,...,p^,q)]l, .(6.17) 
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Figure 3: Illustration of Eq. ( 6.18| ). The left hand side represents H^,{x;pi, . . . ,Pm) expressed in 



terms of V, with the giuon line contracted with q indicated by the curly line with an arrowhead. The 
sum is originally over graphs Gi £ C^, but this is the negative of the sum over graphs Gi € C'_, in 
which the gluon connects to line «, as depicted in the right hand side of the first equality. The remaining 
Green function is W. The insertion then cancels the adjoining propagator. After summing over 
graphs, we are left with a color matrix times the complete m particle tree amplitude. 

Note that we have written the momentum of parton m + 1 as g instead of xpi. Then we take 
the limit as g ^ xpi. This avoids an ambiguous expression of the form 0/0 in the subsequent 
calculation. 

We now note that gauge invariance says something about H. We have inserted a gluon 
line carrying momentum q almost everywhere into tree graphs with m legs. The gluon has 
polarization q^. Gauge invariance tells us that if we inserted this gluon everywhere, that is 
if we summed over the entire set of graphs C", we would get zero. However we have summed 
only over the graphs in C^, leaving out the graphs in C'_. Thus H equals the negative of the 
sum over graphs in C'_. This enables us to calculate H as follows: 

Ha'{x;pi, . . . ,Pm) 

= ~ H.™ X] (l-2;)"^/2g^'u((l-x)pi)[V''(Gi,pi,...,(l-a;)pi,...,pm,g)]^, 

G,&C'_ 

[{1 - x)p^ + qy + lO 
X [W(G-,pi, . . . , (1 - x)pi +q,... ,Pm)]^„ 

= J2 SsfJ'^T^'a"HPi) [W(G-,Jli, . . . . . . ,Pm)]^„ 

= gs^l'T^,^nM{pu...,p^,...,PraT" • (6.18) 

Following the second equals sign, we have displayed the vertex at which the gluon couples to 
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the quark Une i, together with the adjacent propagator. Everything else we call W. Note 
that W is the Green function for an amputated tree graph G'^ with m external partons. The 
partons have flavors {/i, . . . , fm}, just as with our original loop graphs, and they have the 
same momenta as the original partons except that parton i carries momentum (1 — x)pi + q. 
Let us denote the set of all such graphs Cq. The sum over d £ C'_ implies that G[ runs over 
all of Cq. In the next step, we replace ^ by ((1 — x)p^ + (f) next to the spinor. This gives a 
factor ((1 — x)pi + g)^/ ((1 — x)pi + g)^ = 1. After this cancellation, it is safe to take the limit 
q — > xpi. Also, we replace n((l — x)pi) by ^Jl — x u{pi). In the last step, we recognize that 
we have the complete three level amplitude M for the m external particles. This calculation 
is illustrated in Fig. ^. 

We can now insert this result into Eq. ( |6.16D and perform the loop integral to obtain 

CiiG;pi,...,p^r 

G,ieic{G) 

d% fuv{lj,lj -Pi) f^dx f h \ ..err^a 



— dlx-^ 2igsfi'T- 



(27r)M^| + iO)((pi-/i)2 + iO) 7o X V P^■ni. 

X i [(1 - X) - 1] gsH'T^fa" M{pi,...,Pi,..., PmY 



aa' 



= ^^^fll^ ("7 + ^(P^^ ■■■,P^,■■■ ,Pmr ■ (6.19) 

Thus, when we sum over all graphs the collinear subtraction terms associated with an external 
quark line with label i, we get a simple singular factor times the tree level amplitude M. 

For an external antiquark line we get the same result from essentially the same calculation. 
We just have to exchange u{pi) for v{pi), reverse the order of matrix multiplications in Dirac 
spinor space, and replace ^ — > —fi in the Dirac propagator numerators. 

6.2 Gluon line 

In this subsection we calculate the sum over graphs of the collinear subtractions |Ci(G)) in 
the case that the external line i is a gluon line. There are two possibilities. Either the external 
gluon connects to a virtual quark and antiquark, or else it connects to two virtual gluons. 

The case in which the external gluon connects to a virtual quark and antiquark is simple. 
In this case the collinear limit of the numerator is proportional to 

x(l - x)p.^{pi)p. = x{l - x) [-pj^{p^) + 2f.pi-e{p,)\ = , (6.20) 

where £^{p) is the polarization vector of the external gluon. The first term is zero because of 
the pf = 0. The second term is also zero because we consider only physical polarizations of 
the external gluon, so that pi ■ £{pi) = 0. Thus |Ci(G;pi, . . . ,Pm)) is zero when the external 
gluon connects to a virtual quark and antiquark. 
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Figure 4: Illustration of Eq. (6.23). The left hand side represents f^''^{G;x;pi 



defined in Eq. (6^) as the coUinear limit of Z| {Ij—piY times the integrand G{G; h, . 



■ ,Pm)°', which is 

J J Pi J ■ ■ ■ I Pm ) ■ 

Here the integrand Q consists of a polarization vector s'^{pi), a vertex and two propagators times an 
amputated Green function V, as in Eq. ( 6.21 ). As in Fig. |2[ the illustration represents the integrand: 
an integration over the loop momentum is not implied. On the right hand side, there are two terms. 
In each, one gluon line of V is contracted with its momentum vector (xpi) or (1 — x)pi respectively, 
indicated by the arrowhead. The other gluon line is contracted with the original polarization vector 
e''{pi), which is written as e^{{l — x)pi) in the first term and s^{xpi) in the first term. The remaining 
factors in the right hand side of Eq. ( |6.23| ) are indicated. 



We now turn to the case in which the external gluon connects to two virtual gluons, 
g + g ^ g. Our analysis is similar to that of the previous subsection for q + g ^ q- There 
are, however, some subtleties, so we present the argument in some detail. 

By assumption, line i, carrying momentum pi out of the graph, connects to a virtual loop 
in G at a gluon-gluon-gluon vertex. We suppose that the gluon lines in the loop have labels 
j and j + 1. Thus the gluon line with label j carries momentum Ij into the vertex and the 
gluon line with label j + 1 carries momentum —Ij+i = Pi — Ij into the vertex. We write the 
integral for the graph as 

f cfil ~ 

G{G,pi,...,pm)b = J j^;^ g{G;li,. . . ,ln,Pl, - ■ ■ ,Pm)b 

= „^T'^ e''{pi)v!;ll{p,,l,-p,,-l,) 
"'J i2nr (/2 + io)((p,-/^.)2 + iO) 

x[V{Gi,pi,...,pi-lj,...,pm,lj)]'2 ■ (6-21) 

We display a gluon color index b for line i in Q, leaving the other indices on Q unwritten. The 
amplitude V has a color index a corresponding to the gluon j in the loop and a gluon color 
index c corresponding to line j + 1 m the loop. There is an explicit color matrix = ifbac 
connecting the colors of V to the color of Q. The amplitude V also carries vector indices /i, v 
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corresponding to the polarization of the gluon hnes j,j + 1. The polarization vector e'^{pi) 
represents the final state gluon on line i. We have displayed the integration over the loop 
momentum Ij, the propagators for the gluons in the loop and the gluon-gluon-gluon vertex 
V, 

Vaufi{Pi,lj -Pi, -Ij) = gau{lj - '^Pi)^l+ gufiiPi - '2lj)a +g^ia{Pi + lj)u ■ (6.22) 

Everything else is included in the Feynman amplitude V for the amputated tree level graph 
Gi obtained by omitting the propagators j and j ' + 1 and the vertex that attaches these 
propagators to external line i in graph G. Thus the amplitude V is defined by Eq. ( |6.21 ). 
The momentum arguments of V are the external momenta of the original graph, plus the 
momenta Pi — Ij and Ij of the two gluon lines. 

With the notation thus defined, we are ready to calculate. The subtraction term f*-"'^ 
defined in Eq. ( |6.2| ) can be expressed in terms of V as 

ff'^{G;x;pi,.. . ,pm)b = lim Z| [pi - Ijf G{G;li, . . . . . ,pm)b 

= lim {-igs ^J.^)T^c£''{pi) Vav^{pi, Ij - Pi, -Ij) [V(Gi,pi, ...,Pi- Ij, . . . ,Pm, lj)]Z 

Ij >X pi 

2 — X ~ 
= ^Esl^^T^c—^ {xpi)^,£,y{{l - x)pi) \V{Gi,pi, . . . , (1 - x)pi, . . . ,pm,xpi)]'^'^ 

1 + X ~ 
+igs^%aYi:^ ii^ - x)Pi)'^^f^(.^Pi) [^(.Gi, Pi, ... ,{1 - X)pi, ... ,Pm,Xpi)]'^'^ . (6.23) 

There are two terms here. In the first, we contract the polarization index u of the gluon line 
with momentum (1 — x)pi into the index i/ of the polarization vector e^{pi) = £u{{^ — x)pi), 
while we contract the polarization index ji of the gluon line with momentum xpi with the 
index of its momentum vector (xpi)^. In the second term the situation is just reversed, so 
that gluon line with index v and momentum {l — x)pi is contracted with its momentum vector. 
It is useful to represent these two terms diagrammatically as different graphs, depicting the 
gluon line contracted with the polarization vector as a curly line and the gluon line contracted 
with its momentum vector as curly line with an arrowhead, as illustrated in Fig. We can 
denote these two graphs as Gj and G^. 

In the case of the second term, we revise the notation, exchanging fj, ^ u, a ^ c, 
X ^ (1 — x) and exchanging the positions of the xpi and (1 — x)pi arguments of V. Then 
the two terms have the same form except that in one term we have the amplitude V(Gj, . . . ) 
and in the other we have V(G?, . . . ) 

As in the previous subsection, we would now like to sum over graphs. We consider one- 
loop amputated graphs G with m external particles having fiavors {/i, . . . , fm} and momenta 
{pi, . . . ,Pm}- We are considering the collinear subtraction for parton i, which we assume 
here is a gluon. We consider only graphs that have a collinear divergence for external leg i: 
i S Ic{G). We call this class of graphs C. 

The corresponding graphs G\ and Gf are amputated tree graphs with m + 1 external 
legs. The fiavors of the first m external particles are {/i, . . . , fm}, the same as for the graph 
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G. The momenta of these particles are the same as for G except for particle i, which carries 
momentum (1 — x)pi instead of pi. The external particle with index m + 1 is a gluon with 
momentum xpi. Let us call the set of all such graphs C . When we sum over all graphs 
G G C, the graphs Gj together with the graphs Gf include all graphs in C' except for the 
graphs in which gluon m + 1 couples directly to gluon i. These graphs are absent because 
graphs with a self-energy insertion on external line i are not included in C. Let us call the 
set of (m + l)-particle graphs that we do get and the graphs that we don't get C'_. 

The collinear subtraction defined in Eq. ( |6.12| ), summed over graphs G £ C, is 



^ Ci{G;pi,. . . ,pm)b 

f d% fuvilj,lj-P^) f\rs(r ^^'"^'^ V f^(r-r-m 

J (/|+io)((p.-/,)^+io) h h 

J {27ry (/|+iO)((p,-/,.)^ + 10)70 V P^-nJ'^'^''"^ 



...,Pm)-- HU0,PU ....Pm)- ' • • • '^-^ ] ,(6.24) 



x[l — X) X [1 — x) 



where we define 



Hac(x;pi, . . . ,pm) = 

(1-x) lim q^eu{{l-x)pi)[\'{Gi,pi,...,{l-x)pi,...,pm,q)]'t'^ ■ (6.25) 



We now use gauge invariance. We have inserted a gluon line carrying momentum q with 
polarization almost everywhere into tree graphs with m legs. Gauge invariance tells us 
that if we inserted this gluon everywhere, that is if we summed over the entire set of graphs 
C", we would get zero. However we have summed only over the graphs in C^, leaving out the 
graphs in C'_. Thus K equals the negative of the sum over graphs in C'_. This enables us to 
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Figure 5: Illustration of Eq. (6.26). The left hand side represents Hac{x;pi, . . . ,pm) expressed in 
terms of V, with the gluon line contracted with qi indicated by the curly line with an arrowhead. 
The sum is originally over graphs Gi G , but this is the negative of the sum over graphs Gi E C'_ , 
in which the gluon connects to line i, as depicted in the right hand side of the first equality. The 
remaining Green function is W. The insertion then gives two terms, not depicted in the figure. 
One of the terms cancels when summed over graphs, while in the other the propagator adjoining the q^ 
insertion is cancelled. After summing over graphs, we are left with a color matrix times the complete 
m particle tree amplitude. 



calculate H as follows: 
Hac{x;pi, . . . ,Pm) 

= -{l-x) lim q^Su{{l-x)pi)[Y{Gi,pi,...,{l-x)pi,...,pm,q)]'!'^. 

- n ^\ \^ a n'T^ g;.£^((l " ^^^^'^^^((l - x)pi, -(1 - x)p, - Q, q) 

— —[L — X lim > BsfJ- J^cd 77^ ^ ; NO , -n 

G.GCo 

X [W(G-,pi, . . . , (1 - x)pi + q,... ,Pni)]^^a 
= (l-x)^lim |g.M%de.((l-x)K) [WiG[,pi,...,{l-x)p^+q,...,Pm)Ya 

r ;/T" ^ ' '^^^ ~ ^^^'^ ~ ^^^^ ^ Wrr^ n. (I r)v+a 
-gs^i J-cd ((1 _ 3.)^. + g)2 _^ io W(_Gi,]5i,...,(i-2;jp, +g, ...,:PmjJ^j> 

= {l-x)gs^l'T^deu{p^) [W(G^,pi,...,K,...,p^)]^ 

= {l-x) gs^Mipi, . . . ,pi, . . . ,pm)'' . (6.26) 

After the second equals sign, we have displayed the vertex at which the gluon with momentum 
q couples to the gluon line i, together with the adjacent propagator. This multiplies the Green 
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function W for an amputated tree graph G'^ with m external partons. The partons have the 
same flavors as in our original loop graphs, and they have the same momenta except that 
parton i carries momentum (1 — x)pi + q. We denote the set of all such graphs Cq. The sum 
over Gi G G'_ implies that G[ runs over all of Cq. The vector and color indices displayed for 
W are those of the line i. Now, when we evaluate (7^ contracted with the vertex function 
and use (1 — x)pi • — x)pi) = as well as = 0, we get two terms, as indicated on the 
right hand side of the third equality. In the second term, the amplitude W is contracted 
with the momentum carried by the gluon line. After summing over graphs G'^, this term 
vanishes. Finally, we take the limit q xpi and recognize that we have the complete three 
level amplitude Ai for m external particles. This calculation is illustrated in Fig. ^. 

We can now insert this result into Eq. ( |6.24 ) and perform the loop integral to obtain 



^ Ci{G;pi, . . . ,pmf' 



y (27r)'^ (/2+i0)((p,- /,)2+ 10)70 V Pi-nJ 



x) 2 



X 



iPil ■ ■ ■ iPn 



■2^^,2ef d% fuv{ljJj-Pi) 

-'^^^^^ J Wf ii] + io)(b. - h? + io)-^^^^' • • • • • 

^CaJ^^ (-- + Oie)) Mipu . . . . . . ,Pn.)' . (6.27) 



Att r(l - e) V e 

Thus, when we sum over all graphs the collinear subtraction terms associated with an external 
gluon line with label i, we get a simple singular factor times the tree level amplitude M. 

7. Final formulas 

There are two steps to the algorithm that we have outlined here for generating the expressions 
to be used in a numerical calculation of the one-loop graphs for a QCD amplitude. 

The first step applies graph by graph, generating subtractions for each graph. If the graph 
is ultraviolet divergent, one generates a subtraction term for the integrand of the graph. The 
subtraction term matches the integrand when the loop momentum is much larger than the 
external momenta. In this way, the ultraviolet divergence of the integral is removed. The 
subtraction terms are designed to have the same effect as standard MS renormalization. If the 
graph has soft or collinear divergences, or both, then one generates corresponding subtraction 
terms. The subtraction terms match the integrand in the soft and collinear limits of the loop 
momentum. In this way, the soft and collinear divergences of the integral are removed. 

The second step generates a next-to-leading order contribution to the m-parton ampli- 
tude that is proportional to the tree level matrix element (summed over tree graphs). This 
contribution has the form 

l''{e)\Mtree{pi,---,Pm)) , (7.1) 
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where I^(e) is a singular function of the dimensional regularization parameter e and is a linear 
operator on the color space of the amplitude j^Wtree {Pi, ■ ■ ■ iPm) )■ There are three parts to 
I^(e). The first comes from the factors -^/rj for each external line that relate the scattering 
matrix to the amputated Green function. These contributions are given in Eq. ( 3.19| ). There 
is one term for each external line i. The second contribution comes from adding back the 
collinear subtractions. Again there is one contribution for each external line, as given in 
Eq. ( |6.19 ) for quarks and Eq. ( |6.27| ) for gluons. Finally, there is a third contribution that 
comes from adding back the soft subtractions. There is one term for each pair of external 
lines, as given in Eq. ( 5.15| ). The net result is 



(4^) 
An r(l - 



-2pi -pj 



^ m 



7i 



(7.2) 



1=1 



where the ji factors are 



Ig = y - -TRHf 



(7.3) 



The and 1/e poles are all of infrared origin. 

We now turn to the complete calculation of an infrared safe cross section a, which we 
write in the form (following as much as possible the notation of ||l5|) 



^ _ ^LO{m-nj} _|_ ^NLO{m+l-nj} _|_ ^NLO{m-nj} 



(7.4) 



There are a lowest order term and two next-to-leading order terms. We suppose that there 
are n/ initial state particles. (In applications, ni is 0, 1, or 2.) In the LO term there are m 
total particles, rij in the initial state and the remaining m — n/ in the final state. We denote 
this contribution by fjLOjm-n/}^ |^^^ discuss its calculation. The NLO term from real 

parton emission has m + 1 — nj final state particles. This is calculated from the m + 1 — nj 
particle matrix element minus subtractions. One can use, for example, the subtraction scheme 
proposed by Catani and Seymour [^]. We denote this contribution by o-NLO{m+i-n/}^ 
not further discuss its calculation. The remaining NLO term has m — nj final state particles. 
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It is calculated as follows: 



^NLOW-n,} ^ ^^™-n, ^.j-".) (^^ , . . . , |. 



X 



(2vr)4 

2Re|(A1tree (pi, • • • ,Pm) \§{G; {l};pi, ■ ■ ■ ,Pm)) 

-{Mtree {Pl, ■ ■ ■ ,Pm) \TZ{G ] {1} ] pi , . . . ,Pm)) 

X] {Mtree {pi,---,Pm)\Sij{G;{l};pi,...,Pm)) 



J2 (-^trec (pi, ■ ■ ■ ,Pm) \Ci{G] {l};pi, . . . 



ie/c(G) 



x21im(7Wtree(pi,---,Pm)|l^(e)+I^(e)|A^tree(pi,---,Pm)) • (7.5) 



In the first term here, there is a sum over amputated one-loop graphs G. We integrate 
over the phase space d^ for m — nj final state particles with momenta Pm+i, ■ ■ -Prm where 
the initial state particles have momenta pi, . . . ,pni- Next, we supply a measurement Fj 
appropriate to the infrared safe observable that we wish to calculate. ("J" is for "jet".) Then 
there is an integration over the loop momentum I as described in Sec. |3.2| . Now there follow 
four terms inside the loop integration. The first is the integrand G for the graph G times 
the complex conjugate of the tree amplitude for m — nj final state particles. The second is 
{M tree times the renormalization subtraction TZ, which is zero if the virtual loop graph is 
ultraviolet convergent. The renormalization subtraction is defined in Sec. |^. The third term 
contains the soft gluon subtractions S for graph G, one subtraction term for each pair {i, j} of 
external lines that corresponds to a soft divergence. The soft subtraction is defined in Sec. ^. 
The fourth term contains the collinear subtractions C for graph G, one subtraction term for 
each external line i that corresponds to a collinear divergence. The collinear subtraction is 
defined in Sec. |^. 

There is one final term, which has no integral over a loop momentum and involves only 
the tree amplitude. Here we have the singular function I^(e) from Eq. (7.2) that adds back the 
soft and collinear subtractions from the virtual graphs. We also have a singular function I^(e) 
that adds back the soft and collinear subtractions from the real graphs. Naturally, exactly 
what function goes here depends on what subtraction scheme one uses for the real graphs. 
However, the 1/e and poles, whose structure follows from the structure of QCD, cancel 
between I^(e) and I'^(e). (For instance, one can easily check that this cancellation works for 
the scheme of [|l5|, where 1(e) corresponds to our 2I^(e).) Generally a finite contribution 
remains in the limit e — > 0. 
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A. UV counterterms 

In this appendix we give the ultraviolet counterterms for the one-loop propagator and vertex 
corrections in the massless case. We follow the general prescription given in Sec. ^, which 
works in every case except for the gluon self-energy, where we need an extra subtraction that 
was not discussed in Sec. Q in order to handle the quadratic divergence. 

Since the divergent graphs and the corresponding counterterms depend on the form of 
the Feynman rules it is useful to define them. In Feynman gauge we have the following rules: 



Gluon propagator: 



-^D^^ik) = --^5^, , (A.l) 



A;2 + iO 

where a and b are the color indices and k is the momentum carried by this line. 



Quark propagator: 



where a, (3 are the color indices and p is the momentum carried by this line. 

• Ghost propagator: 

where o, h are the color indices and k is the momentum carried by this line. 

• Gluon-quark vertex: 

r^„ = ig,/.Vt^, , (A.4) 

where fj^^ is the fundamental representation of the generators of the gauge group and 
a, a and /3 are the color indices of the gluon, antiquark and quark respectively. 

• Three-gluon vertex: 

G™(pi,P2,P3) = ^■^gs^^'F:l,^V^^^^-^^{p^,p2,p^) , (A.5) 

where F^^^^ = —ifaia2a-i is the adjoint representation of the generators of the gauge 
group and the kinetic part of the three-gluon vertex is 

y'^™(pi,P2,P3) =5^^''Hpi -PzY' +9^'^'{Pz-P2Y' +g^'^'{p2-piY' , (A.6) 
where the pi momenta are outgoing. 
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• Four-gluon vertex: 



010203114 '■Sjsr' V 01 02 0304 \y y ij y ) 

_|_7?f' 7?'' ('„MlM2 „/i3/i4 _ „/il/13„A'2^^4^ 

0104-^ 0203 \y y y y J 



where the Oj are color indices. 



• Gluon-ghost vertex: 

GL(fc) = '^^sfi'FbX , (A.8) 



where a, b, and c are the color indices of the gluon, ghost, and antighost legs respectively 
and k is the momentum of the outgoing ghost. 

We use the standard notation for color factors, Ca = { = 3 ioi SU(3)), Cp = (N^ — 
1)/{2N,), Tr = 1/2. 



Quark self-energy 

b 



l-h 




The quark self-energy graph is the simplest one-loop Green 
function. It has logarithmic and linear UV divergences and we 
can easily define the UV counterterm by applying the prescrip- 
tion introduced in Sec. ^. 

This graph is shown in Fig. ^ and the renormalized quark 
self-energy is defined by the following integral: 



Figure 6: Quark self-energy 
graph. 



dH 
{2ttY 



T.^p{l,p)-Ty^J{l,p) 



where the integrand of the self-energy graph is given by 



af3 



i + \pf + m [{i-\pf + iQ 



and the corresponding ultra-violet counterterm is 



^al3 {hP) - -Igs/^ (-^F0al3- 



{P - //2e-i + iO)^ 
Here e = 2.71828 ... is the base of natural logarithms. 



(A.c 



(A.IO) 



(A.ll) 
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Gluon self-energy 

The one-loop gluon self-energy is sum of the four contributions shown in Fig. (|^ plus an 
ultraviolet subtraction, 



\iiZ{hp)--^Kluvi^^p) 



(A.12) 



Applying the Feynman rules, one finds that the integrand H'^'^ is 



2! (/f + iO)(/i + iO) 

Tr [rhrk] 



(/f + iO)(/i + iO) ' 



(A.13) 



where a, b are the color indices, is the kinematical part of the three-gluon vertex, and 

the loop integral is parametrized by 



1 



l'i=l'- l,P 



ll^ = l^ + ip^ 



(A.14) 



The integrand of the gluon self-energy has quadratic and logarithmic ultraviolet diver- 
gences. As pointed out in Sec. ^ the quadratic divergence requires a subtraction defined by 
a somewhat different prescription than the simple prescription described in general terms in 




Figure 8: One-loop contributions to the quark-gluon vertex. 



Sec. Q and used elsewhere in this appendix. We choose the counterterm to be 

{{d - 2)Ca - iTRUf) (2m^ + 25^ (5^ V - P^P") 



Kluvihp) 



ab 



V)^ + iO) ((/ + ip)' + iO 



3d - 2 Ca(£/''V - P^P"") 2d -A TRnf{g^"'p'^ - p'^p" 



2d - 2 (/2 _ ;^2gl/15 + i0)2 d-l (/2 _ ^2g-l/3 + i0)2 



. (A.15) 



The first two terms eliminate all the quadratic divergences and the integral of them is exactly 
zero in d = 4 — 2e dimensions. The last term cancels the remaining logarithmic divergences. 
The integrand [n^^(/,p) — n^^^y(/,p)] has the familiar tensor structure g'^'^p^ — p^p'^. 

Quark-gluon vertex 

At one-loop level there are two graphs that contribute to the quark-gluon vertex, as shown 
in Fig. p. These graphs have only logarithmic ultraviolet divergences. The corresponding 
renormalized quark-gluon vertex is 



r;i|a(Pl,P2,P3) = igs/^'t^a7M 



Pa 



d^^l 
{2'kY 

dH 
{2t:Y 



Q^.{l,Pi,P2,Pz)-Q''/{i) 
G^{l,Pi,P2,P3)-G'^''{l) 



(A.16) 



The integrands Qn and Gn are derived directly from the Feynman rules, 



Qfi{UPl,P2,P2.) = -igsM^' ( Cf 



Ca 



2 J (z2 + io)(/2 + io)(/i+iO) 



/, N .2 2eCA Ykl^Vf^uaiPl, 13,-12) 

G^,il,Pl,P2,P3) = -Igs/^ 



2 (/2 + io)(/2 + io)(;2 + io) ' 



(A.17) 
(A.18) 
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Figure 9: One-loop corrections to the three-gluon vertex. The {1,2, 3}' denotes the cychc permutation 
of the indices (1,2,3). 



where V^y^ is the kinematical part of the triple gluon vertex and the loop integral was 
parametrized by 



1^ = 1^" + 



P3-P2 



lii = 1'' + 



3 ^ 3 

The corresponding ultraviolet counterterms are 



Pi -P3 l^" = lf' + F2 - Pi 



Q7(0 = -ig.V^(rf-2)(Ci.-^) 
G^^(0 = -ig.V^(d-2) 



Ca\ - 2//, 



2 J {P- /i2e-2 + i0)3 ' 
Ca /27^(l + e) + 2//^ 



2 (/2 _ ^2g-2/3 + i0)3 ■ 



(A.19) 

(A.20) 
(A.21) 



The quark-photon vertex correction can be gotten from the quark-gluon vertex with the 
substitutions Ca 0, Cp ^ I, g ^ —Q\e\ and — > 1, where Q is the quark charge in units 
of |e|. Thus the renormalized quark-photon vertex is 



^f^{Pi,P2,P3) = -iQ|e|7/. - i<3|e| 



Three-gluon vertex 



(27r)« 



Q^il,Pl,P2,P3)-Q]i''{l) 



Ca=0 
Cf=1 



(A.22) 



The three-gluon vertex is sum of the five triangle graphs and three bubble graphs shown in 
Fig. ^ It is useful to group the triangle graphs with a ghost loop together with the triangle 
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graphs with a gluon loop. The other graphs are treated separately. Then the renormalized 
triple gluon vertex is 

+ ig,/x2^F«i„3 , (A.23) 

{1,2,3}' ^ 

where T/^i/^z^is jg h^^q contribution of all the triangle graphs, ^'^I'^^/^a represents one of the 
bubble graphs, and {1,2,3}' denotes the cyclic permutations of the indices (1,2,3). The 
integrand of the triangle graphs is given by 

[i,PuP2,P3)- ig./^ I 2 (Zf + iO)(Zi + iO)(Zi + iO) 

r< . ll^l 7^3 I /Ml 7M2 ]H3 

i^A '3 'i i'2 "T '2 '3 n 



+ 



2 (zf + io) (q + io) iq + io) 



+ ^^"^ (/f + io)(/i + io)(/i + io) j • (^•'"^ 

The loop integral here is parametrized as 

f^^ln^Pj^ , l^ = lf^ + PL_A , i^ = if^ + Pl^ . (A.25) 
With these choices, the ultraviolet counterterm is given by 

where the numerator functions are 

{1,2,3}' ^ 

+ ^^'^^'^^ (7Z2(P2-Pir +2Z-(p2 -Pir^)) (A.27) 

and 

^WM2;^3(;) = /'_4/2^Mi;^2^;^3 + ^^Mi;M2;/^3_^;w;M2(p2-i?ir3 

{1,2,3}' ^ 

+ I9^"''{l^(p2-Pir'-1-{P2-Pi)l''')^ . (A.28) 
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The bubble integral is simpler. In this case the integrand has only a logarithmic diver- 
gence. The integrand is given by 



B>'^>''>''{h,l2,Pi) = iggM'^7<^A " ,,o , ."^ , (A.29) 



9 gfJ-l/J-ipl^'^ _ gf^lt^2pl^i 

4^^ (Zf + iO)(/i + iO) 



with the following parametrization for the loop integral, 



The corresponding counterterm for the bubble integrals is 

4^ (Z2-^2+i0)2 



Four-gluon vertex 

The one-loop correction to the four-gluon vertex is sum of box, triangle, and bubble graphs. 
In the case of the box graphs, it is useful to group the graphs with a ghost loop together with 
the graphs with a gluon loop. We write the renormalized vertex as 

^aia2a3aiR\Pi^P2^P3^P4) — Igs/^ ^ aia2a3a4 



/ (£p [«r(^Pi-P2,P3,P4) > (A.32) 

where the Pi, P2 and P3 are certain sets of permutations of the indices (1,2,3,4), 

Pi = {(1,2,3,4), (1,2,4,3), (1,3,2,4)} , 

P2 = {(1,2,3,4), (1,4,3,2), (1,3,2,4), (4,2,3,1), (3,2,1,4), (4,3,2,1)} , 

P3 = {(1,2,3,4), (1,3,2,4), (1,4,3,2)} . (A.33) 

The integrand of the box diagrams can be obtained by applying the Feynman rules. 



(/f + iO) (;i + io) (/| + iO) (/| + iO) 
. 2 2e 2n/ Tr {t'^H'^H'^H''^) iVf ^^^^^^(?i, /2, ^3, ^4) , . 
^ '^'^ {ll + iO) {il + iO) {il + iO) {il + iO) ' ^ ^ ^ 

where the loop integral is parametrized by 

f = '^P3+P4-P2 f = 1^^- '^P3+P2-P4 

I ju I M+k-p'4 ' Ij, 2pt+k-p'2 ' ^^-^^^ 
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Figure 10: One-loop level correction to the four-point gluon vertex. The Pi, P2, P3 denote the 



permutations of the indices (1,2,3,4) and they are given in Eq. (|A.33| ) 



and the numerator functions Ng and Nq are given by 



^MiM2A.3/.4 (/^ ^ h, h) = V^'^'sih - k, k, -hW'iik -ll,h, -k) 

^v^'^ih - I2, 12, -hW^'^ih - k, h, -h 

_ /Ml /M3 /M4 _ /Ml /M2 /M3 /M4 
M 'l ''2 'S 'S M 'l ''2 ' 

^MiM2M3/.4(^^^ /s, k) = TV [7^V47''Vi7''V27'''/s] • 



(A.36) 
(A.37) 



The integral of the box graphs has only a logarithmic divergence. The ultraviolet counterterm 



- 36 - 



IS 

(/2_^2g-4/3 + i0)4 

+ ig'^ —77Z^ ^4 • A.38 

The integrand of triangle graphs is 

la,a,a„, il,Pl,P2,P3.Pi} - ^S.f (i? + iO) (!| + iO) {i| + iO) ' 

where the loop integral is parametrized by 

li^ = lf^ + P±^PA. , l^ = l^^ + Pl±Pl^P± , ;^ = //^ + ^3 -Pi -P2 _ (A.40) 
3 3 3 

With this parametrization the UV counterterm to the triangle graph is 

^UVa,a,asaAn - Igs/^ _ ^2^-2/9 + jq) ^ " ^^"^^^ 

Applying the Feynman rules, the integrand of the bubble graphs can be found easily, 

1 n n wA'iA'2°/5w/^'<^V3M4 

^aia2a3a4 \hPl-,P2,P3,P4) ^SsP 2\ (|/2 _|_ -q^ ^^2 _|_ -q^ ' l,^."*^; 

where 

With this, the UV counterterm to a bubble graph is 



UVaia^aaa^y-^ -osr^ 2! (^2 _ ^2gl/3 _^ ^q)^ 
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